Higgs-Dilaton cosmology: Are there extra relativistic species? 
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Recent analyses of cosmological data suggest the presence of an extra relativistic component beyond the 
Standard Model content. The Higgs-Dilaton cosmological model predicts the existence of a massless particle 
- the dilaton - associated with the spontaneous symmetry breaking of scale invariance and undetectable by 
any accelerator experiment. Its ultrarelativistic character makes it a suitable candidate for contributing 
to the effective number of light degrees of freedom in the Universe. In this Letter we analyze the dilaton 
production at the (p)reheating stage right after inflation and conclude that no extra relativistic degrees of 
freedom beyond those already present in the Standard Model are expected within the simplest Higgs-Dilaton 
scenario. The elusive dilaton remains thus essentially undetectable by any particle physics experiment or 
cosmological observation. 
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1. Introduction 

Cosmology is entering in a precision era where 
the interplay with particle physics is becoming more 
and more important. A noteworthy example is the 
effective number of light degrees of freedom appea- 
ring in the different extensions of the Standard 
Model (SM) . Any extra radiation component in the 
Universe is usually parametrized, independently of 
its statistics, in terms of an effective number of neu- 
trino species, N cff = + AN eS [fj], where Ajjf 
stands for the number of active neutrinos in the 

sm3. 

The strongest constraints on the effective num- 
ber of neutrino species come from Big Bang Nu- 
cleosynthesis (BBN). A non-standard value of N e g 
increases the expansion rate, which results on an 
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1 In the standard cosmological model with three neutrino 

flavors and zero chemical potential we have N^M = 3 at 

BBN, and N s ¥ = 3.046 at CMB. The small excess in the 
off 

last case with respect to the LEP result [2| is due to the 
entropy transfer between neutrino species and the thermal 
bath during electron-positron annihilation [jj. 



enhancement of the primordial helium abundance. 
Assuming zero lepton asymmetry, the number of 
effective degrees of freedom at BBN turns out to 
be N cS = 3.71+8;^ (68% C.L.) 4]. Note that, al- 
though the existence of extra species is somehow 
favored, the obtained value is still compatible with 
the SM prediction within the 95% C.L. 

Some constraints on N e g can be also obtained 
from the analysis of the Cosmic Microwave Back- 
ground (CMB), although the current limits are sig- 
nificantly weaker than those of BBN. The combined 
analysis of WMAP7 results, Hubble constant mea- 
surements and baryon acoustic oscillations [5j pro- 
vides a value iV off = 4.34±g;|| (68% C.L.). Similar 
and complementary results for smaller CMB scales 
have been also reported by the Atacama Cosmo- 
logy Telescope @ and the South Pole Telescope Q • 
It is interesting to notice the dependence of the e- 
ffective number of neutrino species on the priors 
considered in the different Bayesian analysis exist- 
ing in the literature. While in some references the 
SM value, = 3 Js ruled out at 95% C.L. [I,©], 
in others, such as [UJ, it is not. Besides, if the 
helium abundance obtained from CMB measures is 
taken into account, together with the most precise 
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primordial deuterium abundance [111 ], the BBN re- 
sult becomes perfectly consistent with the SM one 
at the 2a level, N cS = 3.22 ± 0.55 [4 The num- 
ber of extra degrees of freedom is therefore an open 
question to be solved by the Planck satellite, which 
is expected to determine N e g with an accuracy of 
~ 0.3 at 2<T [121, breaking thereby the degeneracies 
with nonzero neutrino masses and dynamical dark 
energy 



13] 



In order to account for the apparent radiation ex- 
cess one can consider several possibilities. It could 
be, for instance, the indication of an extra sterile 
neutrino [3, EH , of relic gravitational waves , 
or arise from other exotic possibilities such as a de- 
caying particle [13, EM El] , the interaction between 
dark energy and dark matter [20 Lor the reheat- 
ing of the neutrino thermal bath 2l[. In this Le- 
tter we will consider a different possibility within 
the minimalistic framework of Higgs-Dilaton cos- 
mology 22J, |23|, [24j . This constitutes an extension 
of the Higgs inflation idea [25| , where the Standard 
Model Higgs doublet H is non-minimally coupled to 
gravity. The novel ingredient of Higgs-Dilaton cos- 
mology is the invariance of the action under scale 
transformations. This extra symmetry leads to the 
absence of any dimensional parameters or scalefl 
The simplest phenomenologically viable theory of 
this kind requires the existence of a new scalar sin- 
glet under the SM gauge group 22], the dilaton %, 
non-minimally coupled to gravity. It corresponds 
to the Goldstone boson associated with the spon- 
taneous symmetry breaking of scale invariance and 
it is therefore massless. This property makes it a 
potential candidate for contributing to the effective 
number of relativistic degrees of freedom at BBN 
and recombination. Indeed, this cosmological test 
seems to be the only available probe for determining 
the existence of the dilaton particle. The coupling 
between the dilaton and all the SM fields (apart 
from the Higgs) is forbidden by quantum numbers, 
which, together with the Goldstone boson nature 



2 In particular it forbids the appearance of a cosmological 
constant term in the action. In Higgs-Dilaton cosmology, 
the late dark energy dominated period of the Universe is 
recovered, at the level of the equations of motion, by repla- 
cing General Relativity with Unimodular Gravity. However, 
both the inflationary and preheating stages considered in 
this Letter take place in field space regions where the dark 
energy contribution is completely negligible. We will thus 
omit this point here. The reader is referred to Ref. l23l 
for details about the phenomenological consequences of Uni- 
modular Gravity in the Higgs-Dilaton scenario. 



of this particle, excludes the possibility of a direct 
detection in an accelerator experiment |22j . 

In this Letter we study the (p)reheating stage in 
Higgs-Dilaton cosmology, paying special attention 
to the dilaton production. This Letter is organized 
as follows. In Sections [5] and [3] we review the Higgs- 
Dilaton inflationary model and show that, given the 
hierarchical structure of the non-minimal couplings 
to gravity, the production of SM particles takes 
place, up to some small corrections, as in the sim- 
plest Higgs inflationary scenario 26, 27, 28j. The 
difference between the two models is described in 
Section 21 where we compute the dilaton produc- 
tion, compare it with the total energy density of 
SM particles at the end of the preheating stage and 
determine its contribution to the effective number 
of relativistic degrees of freedom. The conclusions 
are presented in Section [5] 

2. Higgs-Dilaton inflation 

We start by reviewing the Higgs-Dilaton model 



2J, [23J. In the unitary gauge H T = (0, h/y/2), it 



is described by the following Lagrangian density 



1 
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= ^xX'+Zhh')R--(d X y-U(x,h), (1) 



where we have omitted the part of the SM 
Lagrangian not involving the Higgs potential, 
J^Smlwo] ■ The values of the non-minimal cou- 
plings to gravity can be determined from CMB ob- 
servations and turn out to be hi ghly hierarchical 
(£ x ~ lO- 3 ,^ ~ 10 3 - 10 5 ) [H[29|. Th e scale- 
invariant potential U(x, h) is given by 



Hx 4 



(2) 



with A the self-coupling of the Higgs field. The pa- 
rameters a and ft must be properly tuned in order 
to reproduce the correct hierarchy between the elec- 
troweak, Planck and cosmological constant scales. 
In particular, we must require ft <C a <C 1. The 
smallness of all the couplings involving the dilaton 
field gives rise to an approximate shift symmetry 
X — > x+const., which, as described in Ref. [3(|, has 
important consequences for the analysis of quantum 
effects. For the typical energy scales involved in the 
(p)reheating stage we can safely set a = ft = in 
all the following developments. 

We study here the (p)reheating of the universe 
after Higgs-Dilaton inflation. As emphasized in 
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Ref. 23:], particle production is more easily a- 
nalyzed in the so-called Einstein-frame, where the 
Higgs and dilaton fields are minimally coupled to 
gravity. Performing a conformal redefinition of 
the metric, — fl 2 g^, with conformal factor 
il 2 = Mp 2 (£ x x 2 +ihh 2 ), we obtain 



-9 



Ml - 1 - , 
-JLR--K(x,h) 



U( X ,h). (3) 



Here K(x, h) is a non-canonical kinetic term in the 
basis ((ft 1 , (ft 2 ) = (X,h) 



with 



-Mi 



(4) 



(5) 



and U(x, h) = U(x, h)/tt 4 is the Einstein-frame po- 
tential. In order to diagonalize the kinetic term we 
can make use of the conserved Noether's current a- 
ssociated to scale invariance. It can be easily shown, 
via the homogeneous Friedmann and Klein- Gordon 
equations in the slow-roll approximation, that the 
field combination (1 + 6£ x )x 2 + (1 + 6£h)/i 2 is time- 
independent in the absence of any explicit symme- 
try breaking term 23]. This conservation suggests 
a field redefinition to polar variables in the (h, x) 
plane 



Mr 



loe 



(1 + 6£ x )x 2 + (1 + 6&)fc 2 



M 2 



tan 9 = 



1 + 6& h 
1 + 6£ x X 



, (6) 



(7) 



In terms of the new coordinates, the kinetic term 
Q becomes diagonal, although non-canonical, 



K 



1 + 66 



1 



sin 2 9 + <; cos 2 
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M 2 p t tan 2 9 + fi 
+ £, x cos 2 9 (tan 2 9 + q) 2 

where we have defined 



m\ (8) 



potential depends only on the angular variable 
and it is symmetric around 9 = 



U(9) 



XMp 

4^1 V sin 2 9 + s cos 2 i 



(10) 



It can be easily seen that during the (p)reheating 
stage the values of the oscillating field 9 are much 
larger than the \i parameter for a large number of 
oscillations, tan 2 9 ^> p,. This allows us to neglect 
the \x term in Eq. (j5]) and perform an extra field 
redefinition 

p = j~ 1 r , \(j)\ = 6q— tanh -1 — <,cos9 1 . 

a L J 

(11) 

with 



1 + 64 



and a 



(12) 



The variable (ft is periodic and defined in the 
compact interval (ft £ [—<po,<po], where (fto = 
Mp/a tanh -1 — ? ] corresponds to the value of 
the field (ft at the beginnin g of inflation. As hap- 
pened in Higgs inflation [26|, |27j , the absolute value 
in Eq. (|11[) is required for (ft to maintain the sym- 
metry of the initial 9 field around the minimum of 
the potential. In terms of these variables the La- 
grangian takes a very simple form 



Mi 



=26(0) 



-R 
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-(dp) 2 ~l(d(ft) 



U{(ft) , 
(13 



that has been widely studied in the literature 
32j, [33J, |34[. In our case, the coefficient in front of 
the dilaton kinetic term is given by 

e 2b ^ =<;cosh 2 {a K (<ft - H)] , (14) 

with k the inverse of the reduced Planck mass Mp. 
The inflationary potential (fTU|) becomes 



U(ct>) 



XM4p (l-^cosh 2 [ aK (0 o -|^|)])' 



4^(1-0 



(15) 



„ = & and , = II± 6 ^k 



6 



(i + eya 



(9) 



The dilatonic field r is massless, as corresponds to 
the Goldstone boson associated with the sponta- 
neously broken scale symmetry. The inflationary 



3. SM particle production 

As shown in Fig. [TJ the shape of the Higgs- 
Dilaton potential (fT5"|) clearly resembles that of the 
simplest Higgs inflationary scenario. In spite of the 
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proportional to the small parameter £ x . This al- 



Figure 1: Comparison between the Higgs-Dilaton inflatio- 
nary potential (blue continuous line) obtained from Eg. (I15H 
and the corresponding one for the Higgs inflation mode (red 
dotted line). In spite of the slight differences in the upper 
inflationary region, they nicely agree in the lower part, where 
the (p)reheating stage takes place. Here Uq = AM|;/(4£?) 
and k = Afn 1 . 



slight differences, both of them present an expo- 
nentially flat region for large field values and nicely 
agree for small ones. Indeed, the relation between 
them becomes explicit if we approximate Eq. (|15[) 
for small <f>. The potential around the minimum 
behaves, in a good approximation, as the standard 
chaotic potential 



U(0)^±M 2 HD tf+O(^), 



(16) 



where the higher order corrections can be safely 
neglected after a few oscillations. The curvature 
of the potential, M\ D = (1 + 6£ X )M 2 , coincides, 
up to small corrections, with that of Higgs in- 
flation M 2 — „ c f . The same applies to the 
h 

masses of gauge bosons and fermions. We ob- 
tain (m 2 A j) H D = ™A,f (1 + 6£ x ), with m 2 A f ~ 
aa ^ h p \ 4>\ the Einstein-frame gauge boson and 



fermion masses in the Higgs inflation model 2(j|, 27} • 
Here g — gi , <?2 / cos 9 W and V%yf, for A = W,Z 
bosons and fermions /, respectively. We see there- 
fore that, from the point of view of (p)reheating, 
all the relevant physical scales in Higgs and Higgs- 
Dilaton inflation coincide, up to small corrections 



3 It is interesting to notice at this point that this statement 
is only valid for values of the field such that tan 2 9 2> fJ.- 
For very small values of the angular variable <fi we recover 
the standard A</> 4 Higgs potential. As it happens in Higgs 
inflation [26tl27tl . this region turns out to be extremely small, 
being completely irrelevant for the study of the (p)reheating 
stage. 



lows us to apply the main results of Refs. [26l. \27\ 
to the Higgs-Dilaton case. Let us summarize here 
those results. In Higgs inflation, the SM particles 
are produced through the so-called Combined Pre- 
heating mechanism 26j, 127J • Intermediate W and 
Z bosons are created from the oscillations of the 
Higgs at the bottom of the potential ([To]) , when- 
ever there is a violation of the adiabaticity con- 
dition. While there is no restriction on the num- 
ber of created gauge bosons, the direct production 
of SM fermions by this mechanism is severely res- 
tricted by Fermi-Dirac statistics. The SM fermions 
appear as secondary products of the weak bosons 
created in each zero crossing. Once produced, the 
gauge bosons acquire a large effective mass due to 
the increasing expectation value of the Higgs field 
and decay perturbatively into quarks and leptons 
within a semioscillation of the Higgs field. This 
decay rapidly depletes their occupation numbers 
and postpones the development of parametric re- 
sonance. During the first oscillations, the fraction 
of energy into SM particles is still very small com- 
pared with the energy in the oscillating Higgs field. 
A large number of oscillations (t ~ 300 MZ D ) will 
be needed in order to transfer a significant amount 
of energy into the SM bosons and fermions. The de- 
creasing of the Higgs amplitude due to the expan- 
sion of the Universe eventually reduces the decay 
rate and parametric resonance becomes the domi- 
nant effect. At this point, the gauge bosons start 
to build up their occupation numbers via bosonic 
stimulation and reheating occurs within a few osci- 
llations. Soon afterwards, the Universe is filled with 
the remnant Higgs condensate and a non-thermal 
distribution of fermions and bosons, redshifting as 
radiation and matter respectively. From there on 
until thermalization, the evolution of the system 
is highly non-linear and non-perturbative, which 
makes difficult to make a clear statement about 
the subsequent evolution without the use of numer- 
ical lattice simulations [26J, 127J . However, thermal 
equilibrium is expected to be achieved at a rehea- 
ting temperature T ~ (3 — 15) x 10 13 GeV, much 
above the QCD phase transition scale, Tqod ~ 300 
MeV, due to the large SM couplings (261.127 1. 



4. Dilaton production 

In addition to the SM fields, the Higgs-Dilaton 
inflationary scenario incorporates an extra degree 
of freedom, the dilaton field p. The constancy of 
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the classical background component is of course 
guaranteed by the scale invariance current conser- 
vation, but this reasoning does not apply to the 
corresponding quantum excitations. As suggested 
in Ref. [34 |. these modes can be excited in the 
preheating stage after inflation^ through the non- 
canonical kinetic term in the Einstein-frame La- 
grangian (|13[) . which mixes quantum excitations 
and background solutions. Although the pertur- 
bative dilaton production through this mixing is 
expected to be very small, non-perturbative effects 
might play an important role [34j . In this section we 
estimate the energy density residing in the dilaton 
field at the end of the preheating stage. Let us start 
by considering the linearizecH equations of motion 
for dilaton perturbations 5p k in Fourier space 

Sp k + (3H + 2b) Sp k + ^Sp k = , (17) 

where we have ignored metric perturbations and 
taken into account the constancy of the background 
field p during the Higgs oscillations at the end of 
inflation. The function b — b(<fi) plays the role 
of an additional oscillatory damping term for the 
dilaton perturbations and depends on the absolute 
value of the inflaton field <fi, cf. Eq. §hQ. This 
dependence makes it cumbersome the direct appli- 
cation of the techniques presented in Ref. [34j for 
the study of particle creation in models with non- 
canonical kinetic terms. The field redefinitions used 
there would imply delta functions coming from the 
derivatives of the absolute value, which substan- 
tially complicates the analytic and numerical treat- 
ment of the problem. On the other hand, although 
the rephrasing of Eq. (TP?) as a Hill's equation is 
extremely useful for the understanding of the par- 
ticle creation mechanisms, it is not necessary for a 
precise computation in an expanding background, 
as the one needed to estimate iV ff. For this rea- 
son, we will adopt an alternative approach, dealing 
only with non-singular evolution equations in an 
expanding Universe. Let us rewrite Eq. (IT71) as 

which, after a redefinition of time, dr = a~ 3 e~ 2b dt, 
can be recast in the form of a time-dependent har- 



4 Any dilaton production previous to this stage is com- 
pletely diluted by the inflationary expansion. 

5 The coupled Higgs-Dilaton equations further simplify 
since the Higgs fluctuations are not significantly ampli- 
fied |27|1 , and thus can be treated as decoupled equations. 



monic oscillator 

64 + u,l(r)5 Pk = 0, (19) 

with frequency w^(r) = fc 2 a 4 e 4b . Here the prime 
denotes derivative with respect to the new time r. 
Choosing an initial vacuum state with zero particle 
conteni|j , the number of created dilatons is given 

by 

n k + \ = ^- k (\8p' k f +ul\5p k \ 2 ) , (20) 

and its associated energy density, 
1 f°° 

p x = — / dkk 2 uj k n kl (21) 
Zn Jo 

can be computed numerically by solving Eq. (I19[) , 
together with the background evolution equations. 
The resulting energy density must be compared 
with the energy density in SM particles at the end 
of the preheating stage, C = p\/ P%m- This quantity 
can be easily related to the effective number of light 
degrees of freedom N e g. In order to do that, let us 
note that, once produced, the dilaton particles are 
completely decoupled from the SM particles, being 
its energy density only diluted by the expansion of 
the Universe p°ap = Px a f • Here the subscripts '0' 
and 'f ' stand for the end of the preheating stage and 
the BBN epoch respectively. On the other hand, 
the total entropy of SM particles after thermaliza- 
tion remains constant, sqclq — SfOf . Taking into ac- 
count the relation between the entropy density and 
the number of relativistic degrees of freedom g at 
a given temperature, s — ^gT 3 , the previous ex- 
pression can be rewritten as go^o a o = 3/^/ a /- 
combining this expression with the evolution equa- 
tion for the dilaton energy density described above 
and dividing the result by the energy density stored 
in a single neutrino species, p v = ^g v Tf, we get 

AN eS =(^) = 5o (il) 4/3 C-2.85C. (22) 
\PuJt 9u \9oJ 

In the last equality we have made use of the num- 
ber of SM degrees of freedom at the end of inflation 
(go = 106.75) and at BBN (g { = 10.75). Therefore, 
we see that, in order to have a contribution to the 
effective number of relativistic degrees of freedom 



6 This corresponds to the initial vacuum conditions 
<5p fc (0) = l/v / 2^fc and <5p' fc (0) = -iu k 8p k . 
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within the reach of the Planck satellite, roughly 
a 10% of the energy density at the end of infla- 
tion must be converted into dilatons. Neverthe- 
less, the transferred fraction turns out to be signi- 
ficantly smaller. Evaluating the numerical solution 
of Eq. (THR) at the time at which the energy density 
in SM particles roughly equals the initial en ergy 
density of the inflaton field, r(to ~ 3001/^) [27j |. 
we get C ~ 1CP 7 . The precise value of the C pa- 
rameter weakly depends on the ratio £/, / vA, which 
determines the total energy density available at the 
end of inflation 23j , and it is quite insensitive to the 
particular value of the small non-minimal coupling 

Although the non-perturbative creation of dila- 
tons due to the background field <f> turns out to be 
extremely small, one should consider the possibility 
of producing them as secondary products of the 
Higgs particles created at the preheating and ther- 
malization stages. The Lagrangian (|13[) leads to a 
number of perturbative processes, such as the decay 
of Higgs particles into dilatons ((f) —> pp) or Higgs- 
Higgs scatterings (</></> — >• pp)- The corresponding 
decay rate and cross-section are given respectively 
by 



M 3 HD 
1927rAf|, 



and 



E 2 



576ttM4 



(23) 



where Mg D is the effective Higgs mass in the region 
Mp/(; h < < M P (cf. Eq. (JTSJ)) and E is the Higgs 
energy in the center-of-mass frame. Assuming this 
energy to be of the order of the temperature of the 
thermalized SM plasma, Tq, it can be easily seen 
that the contribution of these processes to the C 
parameter is of order 



Ccx 



M P 



(24) 



and therefore much smaller than the non- 
perturbative contribution found above. 

5. Conclusions 



the model turn out to coincide, up to small correc- 
tions, with those of Higgs inflation. Gauge bosons 
and fermions are therefore produced through the 
so-called Combined Preheating mechanism. On the 
other hand, the production of dilatons could take 
place only through the non-canonical kinetic term 
for the dilaton field, fed by the absolute value of the 
Higgs field. We have shown that the number of non- 
perturbatively created particles can be easily eva- 
luated thanks to a particular time redefinition. The 
dilaton energy density computed this way turns out 
to be extremely small, which is translated into an 
effective number of relativistic degrees of freedom 
very close to the SM one. This result is not modified 
by any of the subsequent perturbative processes in- 
volving the Higgs particles produced at the bottom 
of the potential. We thus conclude that, in spite 
of the potential value of BBN and CMB for test- 
ing the existence of the elusive dilaton particle, its 
subdominant production at the (p)reheating stage 
after inflation makes it completely undetectable by 
any particle physics experiment or cosmological ob- 
servation. The only remnant of the dilaton field is 
a dynamical Dark Energy stage with an equation 
of state close, but slightly different, to that of a 
cosmological constant, which leads to a power-like 
expansion of the Universe in the far future 22j. 
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We have considered particle production in a 
scale-invariant extension of Higgs inflation known 
as Higgs-Dilaton inflation. This model predicts 
the existence of an extra massless particle - the 
dilaton - which might contribute to the effective 
number of light degrees of freedom. After recast- 
ing the problem in the appropriate set of varia- 
bles, all the particle masses and energy scales in 
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